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0.  ABSTRACT 

Let  t  be  the  class  of  all  procedures  for  monitoring  the  variance 
of  a  process  at  level  ct  using  control  charts  based  on  statistics  {T(Y^), 

T(Y2),...}  from  independent  random  samples.  Suppose  the  control  variance 

2  2 
Oq  is  known.  Under  Gaussian  assumptions  the  S  chart  using  the  sample 

variance  is  shown  to  be  optimal  in  the  class  r ^  in  that  its  run  length 

is  stochastically  smallest  under  both  stationary  and  drifting  processes 

not  in  control.  Weaker  properties  are  given  in  terms  of  stochastic 
2 

bounds  when  aQ  is  not  known  and  instead  is  estimated  in  a  base  period 
using  the  sample  variance. 


Acknowledgment .  This  research  was  supported  in  part  by  the  U.S.  Army 
Research  Office  under  grant  DAAG-29-78-G-0172. 

Key  words  and  phrases.  Statistical  quality  control,  process  variance, 
2 

S  charts,  run  lengths,  stochastic  ordering. 


ii 


1 .  INTRODUCTION 


Various  control  charts  are  used  in  monitoring  the  parameters  of  a 

production  process.  For  monitoring  the  variance  of  product  quality  the 

2 

standard  procedures  are  the  R  and  S  charts  using  the  sample  range  and 

the  sample  variance,  respectively.  In  these  charts  successive  values 

of  a  statistic  are  plotted  against  time,  and  a  chart  signals  at  level  a 

that  the  monitored  process  is  not  in  control  when  the  statistic  exceeds 

a  control  limit  c  .  The  standard  assumptions  are  that  (i)  the  process 

2  2 

is  stationary  in  its  variance  with  value  a  ,  (ii)  the  control  value 

is  known,  (iii)  independent  random  samples  of  n  observations  each  are 

taken  on  successive  monitoring  occasions,  and  (iv)  the  sampled  process 

2 

is  Gaussian.  For  the  R  and  S  charts  the  statistics  are  {R^/j^  ;  i  =  I , 

2  2  2 
2,...}  and  {S^/aQ ;  i  ■  1,2,...},  respectively,  where  and  are  the 

sample  range  and  the  sample  variance  of  the  ith  sample.  An  account  of 

these  procedures  is  given  in  Duncan  (1974). 

On  each  sampling  occasion  the  problem  is  to  test  at  level  a  the 

2  2  2  2  2 
hypothesis  H:a  *  against  A:a  >  a^,  for  which  S  is  known  to  have 

greater  power  than  R.  The  usual  notions  of  level  and  power  do  not  carry 
over  to  the  entire  monitoring  period,  however,  as  the  standard  charts 
eventually  signal  with  unit  probability  even  when  the  process  is  in 
control.  What  is  required  is  that  a  chart  should  signal  infrequently 
when  the  process  is  in  control  and  more  frequently  otherwise.  Such  prop¬ 
erties  are  summarized  in  the  distribution  of  the  run  length  of  a  chart, 
i.e.  the  number  of  successive  samples  taken  before  the  chart  signals. 

Run  lengths  accordingly  provide  a  mechanism  for  comparing  two  charts, 


for  tf  the  run  length  of  one  chart  is  stochastically  smaller  than  another 

when  the  process  is  not  in  control,  then  that  chart  will  tend  to  signal 

more  quickly.  Under  standard  assumptions  the  run  lengths  of  the  R  and 
2 

S  charts  have  geometric  distributions,  each  with  parameter  equal  to 
the  probability  of  exceeding  the  control  limit  on  any  sampling  occasion; 
otherwise  the  distributions  are  more  complicated. 

2 

Although  tradition  mandates  use  of  the  R  or  S  chart,  clearly  any 

of  numerous  statistics  could  be  used  in  monitoring  the  variance  of  a 

2 

Gaussian  process.  Here  we  show  that  the  known  optimality  of  S  in  test- 

2 

ing  H  against  A  on  any  sampling  occasion  carries  over  to  the  S  chart  in 

the  following  sense.  Among  all  a-level  monitoring  procedures,  the  run 
2 

length  of  the  S  chart  is  stochastically  minimal  when  the  process  is  not 

in  control.  Weaker  results  provide  stochastically  smallest  bounds  when 
2 

Oq  is  not  known  and  instead  is  estimated  using  the  sample  variance  from 
a  base  period. 

2.  THE  MAIN  RESULTS 

Let  {YpY^,...}  be  n-dimensional  observations  taken  on  successive 
occasions;  let  {T(Y^),  TO^),...}  be  the  values  of  a  statistic  used  in 
monitoring  the  variance;  and  let  N^,  be  the  run  length  of  the  correspond¬ 
ing  chart.  Suppose  the  chart  signals  at  level  a  whenever  T(Y  )  t  B, 

•  ~  1 

i.e.  P(T(Y^)  e  B)  *  a  when  the  process  is  in  control.  It  is  clear  that 

P(Nt  >  k)  =  P(T(Yl)  e  8C,  ...,  T(Yk)  e  BC)  (2.1) 

c 

with  B  as  the  complement  of  B.  Under  the  assumption  that  {T(Y^),  T(Y2) , 
...}  are  mutually  independent,  (2.1)  becomes 


3 


k 

P(Nt  >  k)  -  .^(1  -  B±)  (2.2) 

with  3^  ■  P(T(Y^)  s  B).  For  nonstationary  processes  probabilities  of 
the  type  (2.2)  depend  on  a  countable  sequence  3  =  (S^»  Sj*  •••)  which 
we  call  the  drifting  parameters  of  the  process.  Accordingly,  let  F^,(u; 

3)  »  1  -  P(NT  >  u)  be  the  cdf  of  under  drifting.  For  stationary 

processes  (T^),  T^),  ...}  are  distributed  identically,  in  which  i 

Ic 

case  the  cdf  of  N_  is  the  geometric  distribution  G(k;S)  m  1  -  (1  -  8)  j 

with  3  =  P(Nt  £  B) ,  i.e.  F(u;8)  *  G(u;3)  with  3  *  ( S ,  S,  ...).  Let  V 

be  the  collection  of  all  sequences  5  *  (5^,  ...)  such  that  6^  £  (0,1) 

for  all  i,  and  let  {F^,(u;S);  3  e  V}  be  the  family  of  possible  run-length 

distributions  associated  with  a  given  chart.  Any  such  family  is  said 

to  decrease  stochastically  in  8  if,  for  any  3  and  3*  in  V  such  that  8*  >_ 

3.^  for  all  i,  F(u;  3*)  >_  F(u;  8)  for  every  u  >  0.  From  the  form  of 

(2.2)  it  is  clear  that  {Ft(uj  3);  3  c  V)  decreases  stochastically  in  3 

for  any  chart  of  the  type  considered  here. 

Let  t ^  be  the  class  of  all  a-level  procedures  for  monitoring  the 

variance  of  a  process  using  statistics  (T(Y^),  KY^),  ...}.  Our  result 

2 

for  this  class  is  that  under  Gaussian  assumptions  the  S  chart  is  opti¬ 
mal  in  that  its  run  length  is  stochastically  minimal  when  the  process 
is  not  in  control. 

THEOREM  1.  In  the  class  of  all  procedures  for  monitoring  the  variance 

of  a  drifting  Gaussian  process  at  level  a  using  statistics  of  the  type 

2 

(T(Y^),  T^),  the  S  chart  is  optimal  in  that  its  run  length  is 

stochastically  minimal . 


Proof.  Consider  any  procedure  T  in  together  with  its  run-length  dis¬ 


tribution  F^,(u;  8),  and  let  Fg(u;  §*)  be  the  run-length  distribution  of 

2  2  2 
the  S  chart.  Considered  as  procedures  for  testing  H:a  =  Oq  against 

2  2 

A: a  >  Oq  on  each  sampling  occasion,  under  Gaussian  assumptions  the  test 
2  2 

using  is  uniformly  most  powerful  among  procedures  in  (cf.  Leh¬ 

mann  (1959),  page  97),  thus  assuring  that  3*  >_  6^  for  all  i.  Observe 
that  F^,(u;  8)  *  F(u;  8)  and  Fg(u;  8*)  =  F(u;  8*)  have  the  same  form 
apart  from  their  parameters,  and  that  { F ( t ;  3);  8  t  V}  decreases  stochas¬ 
tically  in  8) •  It  follows  that 


Fs(u;  §*)  >  FT(u;  §)  (2.3) 

for  every  u  >  0  and  T  e  t^,  which  is  equivalent  to  the  assertion  of  the 
theorem. 


2 

Somewhat  weaker  results  can  be  shown  for  the  case  that  oQ  is  un- 

2 

known.  A  common  practice  is  to  estimate  Oq  in  a  base  period  when  the 

2 

process  is  in  control  and  to  use  this  estimate  in  lieu  of  Oq  in  main- 
2 

taining  the  R  and  S  charts.  A  similar  approach  may  be  taken  with  other 


charts.  Accordingly,  let  Yq  be  a  vector  of  m  observations  taken  when 

the  process  is  in  control;  define  {Y^,  Y^,  . ..)  as  before;  and  let 

{T(Y^,  Yg) 5  i  *  1,2,...}  be  the  values  of  a  statistic  to  be  used  in  a 

typical  procedure  for  monitoring  the  process  variance.  For  example, 

2  2  2 

T(Y^,  Yq)  takes  the  values  R^/Sq  and  S^/Sq  in  the  modified  R  and  S 

2 

charts  using  the  sample  variance  Sq  from  Yq. 

Suppose  the  typical  chart  signals  at  level  a  whenever  T(Y^,  Yq)  t 
B,  i.e.  P(T(Y^,  Yq)  s  B)  *  ot  when  the  process  is  in  control,  and  denote 


5 


by  the  run  length  of  this  chart.  Corresponding  to  (2.1)  we  have 

P(Nt  >  k)  =  P(T(Yr  YQ)  e  BC,  ....  T(Y^,  YQ)  e  BC) .  (2.4) 

The  run  lengths  for  such  charts  generally  are  not  geometric  even  for 

stationary  processes  owing  to  dependencies  among  {T(Yi#  Yq) ;  i  =  1,2,... 

However,  when  the  process  is  stationary  but  not  necessarily  in  control, 

the  statistics  {T(Y^,  Yq) ,  ....  T(Y^,  Yq)  }  are  distributed  identically. 

c  t  k, 

The  convexity  of  the  conditional  probability  [P(T(Yi>  Y  )  e  B  |Yq)] 
as  a  function  of  Yq  and  Jensen's  inequality  give  the  bound 

P(Nt  >  k)  >  (1  -  S)k  (2.5) 

with  0  =  P(T(YiS  Y  )  e  3)  for  all  i  =  l,2,...,k  under  stationarity , 
which  we  henceforth  assume.  Let  FT(u;  8)  be  the  distribution  of  with 
parameter  3  under  a  stationary  process.  Then  (2.5)  assures  that  the 
geometric  distribution  G(u;  8)  not  only  approximates  FT(u;  3),  but  that 
it  is  a  stochastic  lower  bound  with  G(u;  8)  ^  FT(u;  3)  for  every  u  >  0. 
Moreover,  equality  is  attained  in  the  limit  as  m  -*■  °°,  so  that  (2.5)  is 
the  best  bound  which  holds  for  ail  m. 

Let  t 2  be  the  class  of  all  a-level  procedures  for  monitoring  the 
process  variance  using  statistics  {T(Y^,  Yq) ;  i  =  1,2,...}  such  that 
P(T(Y^,  Yq)  e  B)  >  a  when  a ^  >  Oq.  The  modified  chart  using  S^/Sq 
has  the  property  that  the  geometric  approximation  to  its  run  length  is 
stochastically  minimal  among  such  approximations  for  all  procedures  in 
when  the  process  is  stationary  but  not  in  control.  Although  somewhat 
weaker  than  Theorem  1,  these  results  are  summarized  as  follows. 


»**«**•*- 


THEOREM  2.  Let  be  the  class  of  all  a-level  procedures  for  monitoring 

a  process  variance  using  statistics  (T(Y^,  Yq) ;  i  =  1,2,...}  such  that 

2  2 

P(T(Y^,  Y  )  s  B)  >  ct  when  a  >  aQ.  If  the  process  is  a  stationary  Gaus¬ 
sian  process  not  in  control,  then  the  approximating  geometric  distribu- 
tion  for  the  modified  S  chart  is  stochastically  minimal  among  such 
approximations  for  all  procedures  in  x^. 

Proof.  For  any  procedure  in  let  G(u;  8)  be  the  geometric  bound  for 

the  distribution  of  its  run  length,  and  let  G(u;  6*)  be  the  correspond- 

2  2  2... 
ing  bound  for  the  modified  S  chart  using  TCY^,  Yq)  =  S^/Sq.  Considered 

2  2.  22 

as  normal-theory  procedures  for  testing  H:a  =  aQ  against  A:  a  >  aQ 

on  each  sampling  occasion,  the  class  x^  consists  of  unbiased  tests.  The 
2  2 

test  using  STVSq  is  uniformly  most  powerful  in  this  class  (cf.  Lehmann 
(1959),  page  169),  thus  assuring  that  3*  >.  8-  The  remainder  of  the 
proof  parallels  that  of  Theorem  1. 

2 

In  conclusion.  Theorem  1  establishes  the  strict  optimality  of  S 
charts  in  the  class  x^.  Theorem  2  supplies  a  weakly  optimal  property 
of  modified  S2  charts  which  holds  for  all  m  and  is  strict  in  the  limit 
as  m  «.  The  interpretation  is  that  the  guaranteed  stochastic  lower 
bound  is  minimal  among  such  lower  bounds  for  all  procedures  in  X£- 
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